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This talk summarizes the recent developments in the evaluation of the leading order hadronic contributions to 
the running of the QED fine structure constant a{s), at s = M|, and to the anomalous magnetic moment of the 
muon (p — 2)^. The accuracy of the theoretical prediction of these observables is limited by the uncertainties on the 
hadronic contributions. Significant improvement has been achieved in a series of new analyses which is presented 
historically in three steps: (I), use of r spectral functions in addition to e^e~ cross sections, (II), extended use 
of perturbative QCD and (III), application of QCD sum rule techniques. The most precise values obtained are: 
Aahad(A/|) = (276.3 ± 1.6) x 10"", yielding a-^(M|) = 128.933 ± 0.021, and a^""^ = (692.4 ± 6.2) x 10"^° with 
which one finds for the complete Standard Model prediction = (11 659 159.6 ± 6.7) x lO^^''. For the electron 
{g - 2)e, the hadronic contribution is a^="^ = (187.5 ± 1.8) x 10"". 



1. INTRODUCTION 

The running of the QED fine structure con- 
stant a{s) and the anomalous magnetic moment 
of the muon are observables for which the theo- 
retical precision is limited by second order loop 
effects from hadronic vacuum polarization. Both 
quantities are related via dispersion relations to 
the hadronic production rate in e+e^ annihila- 
tion, 

Ct,ot(e^e^ — > hadrons) 



R{s) 



o■o(e^ 



(1) 



with (To(e"'"e~ — * A^^M ) — 47ra^/(3s). While far 
from quark thresholds and at sufficiently high en- 
ergy y/s, R{s) can be predicted by perturbative 
QCD, theory fails when resonances occur, i.e., lo- 
cal quark-hadron duality is broken. Fortunately, 
one can circumvent this drawback by using e+e^ 
annihilation data for i?(s) and, as proposed in 
Ref. hadronic r decays benefitting from the 
largely conserved vector current (CVC), to re- 
place theory in the critical energy regions. 

There is a strong interest in the electroweak 
phenomenology to reduce the uncertainty on 
a(M|) which used to be a serious limit to progress 
in the determination of the Higgs mass from ra- 
diative corrections in the Standard Model. Ta- 
ble gives the uncertainties of the experimen- 



tal and theoretical input expressed as errors on 
sin^^w Using the former value ^ for q;(M|), 
the dominant uncertainties stem from the exper- 
imental sin^^w determination and from the run- 
ning fine structure constant. Thus, any useful 
experimental amelioration on sin^^w requires a 
better precision of a(M|), i.e., an improved de- 
termination of its hadronic contribution. 

The anomalous magnetic moment a^^ = {g ~ 
2)/2 of the muon is experimentally and theoret- 
ically known to very high accuracy. In addition, 
the contribution of heavier objects to relative 
to the anomalous moment of the electron scales 
as (m^/me)^ ~ 4 x 10"*. These properties allow 
an extremely sensitive test of the validity of the 
electroweak theory. The present value from the 
combined /i"*" and ii~ measurements 



(11 659 230 ±85) x 10 



-10 



(2) 



is expected to be improved to a precision of 
at least 4 x 10^^° by the E821 experiment at 
Brookhaven |^,H). Again, the precision of the 
theoretical prediction of is limited by the con- 
tribution from hadronic vacuum polarization de- 
termined analogously to a(M|) by evaluating a 
dispersion integral using e+e" cross sections and 
perturbative QCD. 
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Table 1 

Uncertainties of the electroweak input expressed 
in terms of Asin^6'w (xlO~'*). Downward arrows 
indicate future experimental improvement. 



is evaluated using the principle value integration 
technique. 

3. MUON MAGNETIC ANOMALY 



Input Asin 6w 



Uncertainty / Source 



Exp. 18 (LEP+SLD) y i 

a(Mz) 23 Aq-1(A/|) = 0.09 11 

mt 15 Am* = 5.0 GeV (CDF+DO)^§ j 

Theory 5-10 2-loop EW prediction 

Mh 150 65 - 1000 GeV 



2. RUNNING OF THE QED FINE 
STRUCTURE CONSTANT 

The running of the electromagnetic fine 
structure constant a{s) is governed by the 
renormalized vacuum polarization function, 
n-,(s). For the spin 1 photon, n.y(s) is 
given by the Fourier transform of the time- 
ordered product of the electromagnetic currents 
j^^{s) in the vacuum {q^'q'' - q'^g'"')U^{q'^) = 
»/d4xe'«-(0|Ta^^(a:)j-^(0))|0). With Aa{s) = 
-47raRe[n.^(s) - n.^(0)] and Aa{s) 
Aaiep(s) -I- Aahad(s), which subdivides the run- 
ning contributions into a leptonic and a hadronic 
part, one has 



a(s) 



.(0) 



1 - Aaicp(s) - Aahad(s) 



(3) 



where 47ra(0) is the square of the electron charge 
in the long- wavelength Thomson limit. 

For the case of interest, s = M|, the leptonic 
contribution at three-loop order has been calcu- 
lated to be [01 



Aaiop(Af|) = 314.97686 x 10" 



(4) 



Using analyticity and unitarity, the dispersion in- 
tegral for the contribution from hadronic vacuum 
polarization reads fictl 



Aahad(M|) 



3n 



s(s - 



Mi) 



le 



,(5) 



and, employing the identity l/{x' — x — ie)e^o — 
P{l/(a;' — x)} + iTrS{x' — x), the above integral 



It is convenient to separate the Standard Model 
prediction for the anomalous magnetic moment 
of the muon, = (g — 2)^/2, into its different 
contributions, 

«r - aQ^°-fa;i^^-fa-^ (6) 

where a^™ = (11658 470.6 ± 0.2) x 10"^° is the 
pure electromagnetic contribution (see Ref. ||l^ 
and references therein), a^^'^ is the contribution 
from hadronic vacuum polarization, and aj^'^^^ = 
(15.1 ± 0.4) X 10^^" accounts for corrections due 
to exchange of the weak interacting bosons up to 
two loops . 

Equivalently to AQ;had(.^^|), by virtue of the 
analyticity of the vacuum polarization correlator, 
the contribution of the hadronic vacuum polar- 
ization to can be calculated via the dispersion 
integral |l|] 



,had 



(7) 



where K{s) denotes the QED kernel |g5| 



K{s) 



x^ [ 1 



X 

~2 



+ (l + x)' 1 + 



1 



X ln(l + x) - a; + 



|i±4x2lnx(3) 

(1 - x) 



with x = (1 - /5^)/(l + /3^) and = (1 - 
4to^/s)^/^. The function K(^s) decreases mono- 
tonically with increasing s. It gives a strong 
weight to the low energy part of the integral (|^). 
About 92% of the total contribution to c^!^'^ is 
accumulated at cm. energies -^/s below 1.8 GeV 
and 72% of a^^^ is covered by the two-pion final 
state which is dominated by the p(770) reso- 
nance. Data from vector hadronic r decays pub- 
lished by the ALEPH Collaboration jl^] provide 
a precise spectrum for the two-pion final state 
as well as new input for the lesser known four- 
pion final states. This new information improves 
significantly the precision of the determina- 
tion [§. 
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4. IMPROVEMENT IN THREE STEPS 

A very detailed and rigorous evaluation of both 
a(M|) and was performed by S. Eidelman 

and F. Jegerlehner in 1995 ||] which since then is 
frequently used as standard reference. In their 
numerical calculation of the integrals (^) and , 
the authors use exclusive e+e~ hadrons cross 
section measurements below 2 GeV cm. en- 
ergy, inclusive R measurements up to 40 GeV and 
finally perturbative QCD above 40 GeV. Their 
results to which I will later refer are 

Aahad(M|) = (279.7 ±6.5) x lO'^ , 

al""^ = (702.4 ± 15.3) x W^'^ . (9) 

Due to improvements on the electroweak ex- 
perimental side these theoretical evaluations are 
insufficient for present needs. Fortunately, new 
data and a better understanding of the underly- 
ing QCD phenomena led to new and significantly 
more accurate determinations of the hadronic 
contributions to both observables. 

(I) Addition of precise r data 

Using the conserved vector current (CVC) it 
was shown in Ref. jl] that the addition of precise 
T spectral functions, in particular of the t~ — > 
TT^n'^h'T- channel, to the e^^e^ annihilation cross 
section measurements improves the low-energy 
evaluation of the integrals (^) and (|^). Hadronic r 
decays into ud' isovector final states occur via ex- 
change of a virtual boson and have therefore 
contributions from vector and axial-vector cur- 
rents. On the contrary, final states produced via 
photon exchange in e^e^ annihilation are always 
vector but have isovector and isoscalar parts. The 
CVC relation between the vector two-pion t spec- 
tral function vj=i{t — > tttt^ and the corre- 
sponding isovector e~^e~ cross section at energy- 
squared s reads 

i=if + - + 47ra^(0) n 

s 

where vj^i(t — > tttt^ Vr) is essentially the 
hadronic invariant mass spectrum normalized to 
the two-pion branching ratio and corrected by 
a kinematic factor appropriate to r decays with 
hadronic spin J = 1 . The two-pion cross sec- 
tions (incl. the r contribution) in different energy 



regions are depicted in Fig. y. Excellent agree- 
ment between r and e+e~ data is observed. For 
the four pion final states, isospin rotations must 
be performed to relate the respective r final states 
to the corresponding e+e~ topologies [p^ . 

Effects from SU{2) violation 

Hadronic spectral functions from t decays are 
directly related to the isovector vacuum polar- 
ization currents when isospin invariance (CVC) 
and unitarity hold. For this purpose one has 
to worry whether the breakdown of CVC due 
to quark mass effects (m„ ^ rrid generat- 
ing dp^J'^ ^ (niu — rrid) for a charge-changing 
hadronic current J'^ between u and d quarks) 
or unknown isospin- violating electromagnetic de- 
cays have non-negligible contributions within the 
present accuracy. Expected deviations from CVC 
due to so-called second class currents as, e.g., 
the decay — )■ tt^ 77 where the corresponding 
e+e^ final state n'^ri (C=-|-l) is strictly forbid- 
den, have estimated branching fractions of the 
order of (m„ - maY /m^ ~ 10^^, while 
the experimental upper limit amounts to B(t — > 
TT-'qvr) < 1.4 X lO--* |l| at 95% CL. SU{2) 
symmetry breaking caused by electromagnetic in- 
teractions can occur in the p^-p^ masses and 
widths. Hadronic contributions to the p^-p^ 
width difference are expected to be much smaller 
since they are proportional to (m^ — rudY . The 
total expected SU{2) violation in the p width 
is estimated in Ref. ^ to be iTp± — rpo)/rp = 
(2.8 ± 3.9) X 10^^, yielding the corrections 

Sal^'^ - -(1.3 ±2.0) X 10"^" , 
JA4l^,(Af|) = -(0.09 ± 0.12) X 10-4 , (10) 

to the respective dispersion integrals when using 
the T~ — > tt^tt^z/t- spectral function in addition to 
e+e" data. 

Evaluation of the dispersion integrals ^) and (Qj 
Details about the non-trivial task of evaluating 
in a coherent way numerical integrals over data 
points which have statistical and correlated sys- 
tematic errors between measurements and exper- 
iments are given in Ref. [Q. The procedure is 
based on an analytical minimization, taking 
into account all initial correlations, and it pro- 
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vides the averages and the covariances of the cross 
sections from different experiments contributing 
to a certain final state in a given range of cm. 
energies. One then apphes the trapezoidal rule 
for the numerical integration of the dispersion 
integrals (||) and (0), i.e., the integration range 
is subdivided into sufficiently small energy steps 
and for each of these steps the corresponding co- 
variances (where additional correlations induced 
by the trapezoidal rule have to be taken into ac- 
count) are calculated. This procedure yields error 
envelopes between adjacent measurements as de- 
picted by the shaded bands in Fig. |^. 



are found to be 0] 

Aahad(M|) = (281.0 ±6.2) x 10"^ , 

ajj^<i = (701.1 ± 9.4) X 10^1° , (11) 

with an improvement for a^^'^ of about 40% 
compared to the previous evaluation (|^), while 
there is only a marginal improvement of Aq;(M|) 
for which the dominant uncertainties stem from 
higher energies. 




Figure 1. Two-pion cross section as a function of the cm. energy-squared. The band represents the 
result of the averaging procedure described in the text within the diagonal errors. The lower left hand 
plot shows the chiral expansion of the two-pion cross section used (see Ref. Q ) . 



Results 

With the inclusion of the r vector spectral func- 
tions, the hadronic contributions to a(i\f|) and to 



(II) Extended theoretical approach 

The above analysis shows that in order to im- 
prove the precision on Aa{AI^), a more accurate 
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determination of the hadronic cross section be- 
tween 2 GeV and 10 GeV is needed. On the 
experimental side there are ongoing R measure- 
ments performed by the BES Collaboration [p^ . 
On the other hand, QCD analyses using hadronic 
r decays performed by ALEPH and re- 

cently by OPAL Q revealed excellent applica- 
bility of the Wilson Operator Product Expansion 
(OPE) m (also called SVZ approach ||), orga- 
nizing perturbative and nonperturbative contri- 
butions to a physical observable using the concept 
of global quark-hadron duality, at the scale of the 
r mass, ~ 1.8 GeV. Using moments of spec- 
tral functions, dimensional nonperturbative oper- 
ators contributing to the r hadronic width have 
been fitted simultaneously and turned out to be 
small. This encouraged the authors of Ref. 
to apply a similar approach based on spectral 
moments to determine the size of the nonpertur- 
bative contributions to integrals over total cross 
sections in e"'"e~ annihilation, and to figure out 
whether or not the OPE, i.e., global duahty is a 
valid approach at relatively low energies. 

Theoretical prediction of R{s) 

The optical theorem relates the total hadronic 
cross section in e~^e~ annihilation, R(so), at a 
given energy-squared, so, to the absorptive part 
of the photon vacuum polarization correlator 



i?(so) = 127rImn(so + ie) 



(12) 



Perturbative QCD predictions up to next-to- 
next-to leading order (a^) as well as second order 
quark mass corrections far from the production 
threshold and the first order dimension D = A 
nonperturbative term are available for the Adler 
ZJ-function [ p5[ , which is the logarithmic deriva- 
tive of the correlator H, carrying all physical in- 
formation: 



D{s) = -127r^s 



2 dll{s) 



ds 



This yields the relation 
27ri / s 



-R(so) 



(13) 



(14) 



where the contour integral runs counter-clockwise 
around the circle from s = sq + ie to s = sq — ie. 



The Adler function is given by 
Df,i-.s)^NcJ2Q}i^ 



f 

1 + do tts + di al + d2 al 

(6 + 28as + (295.1 - 12.3 n /) a^) 



m'j{s) 



+ 



18 



+ Stt^ (1 - a,) 



{mfqfqf) 



327r2 (mkqkqk) 

k 

12.m + 16.^^^ 



(15) 



where additional logarithms occur when ji^ ^ s, 
^ being the renormalization scaleQ and = 
^i^^. The coefficients of the massless per- 
turbative part are do = 1, c'l = 1.9857 — 
0.1153 n/, ^2 = ^2+ /3g7rV48 with /3o = 11 - 
2n//3, d2 = -6.6368 - 1.2001 n/ - 0.0052 - 

1.2395(J2fQf?/NcJ2fQ} and Uf being the 
number of involved quark flavours. The non- 
perturbative operators in Eq. ( [l5| ) are the gluon 
condensate, {(as/'r^)GG), and the quark conden- 
sates, {mfqfqf). The latter obey approximately 
the PCAC relations 



{rriu + m,d){uu + dd) ~ -2/Jto^ , 



ms{ss) 



(16) 



with the pion decay constant /^r = (92.4 ± 
0.26) McV [|l8|. In the_chiral hmit the equations 
/tt = Jk and {uu) = (dd) = (ss) hold. The com- 
plete dimension D — 6 and D = 8 operators are 
parametrized phenomenologically in Eq. (jlj) us- 
ing the saturated vacuum expectation values (Og) 
and (Os), respectively. 

Although the theoretical prediction of R us- 
ing Eqs. (|l^) and ( |l5| ) assumes local duality 

^ The negative energy-squared in D(—s) of Eq. is 
introduced when continuing the Adler function from the 
spacehke EucHdean space, where it is originally defined, 
to the timelike Minkowski space by virtue of its analyticity 
property. 
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and therefore suffers from unpredicted low-energy 
resonance oscillations, the following integration, 
Eqs. (||)/(0), turns the hypothesis into global du- 
ality, i.e., the nonperturbative oscillations are av- 
eraged over the energy spectrum. However, a sys- 
tematic uncertainty is introduced through the cut 
at explicitly 1.8 GeV so that non- vanishing oscil- 
lations could give rise to a bias after integration. 
The associated systematic error is estimated in 
Ref. by means of fitting different oscillating 
curves to the data around the cut region, yielding 
the error estimates A(Aahad(-M|)) = 0.15 x 10""' 
and Aa^'^'^ = 0.24 x 10"^", from the comparison 
of the integral over the oscillating simulated data 
to the OPE prediction. These numbers are added 
as systematic uncertainties to the corresponding 
low-energy integrals. 

In asymptotic energy regions we use the for- 
mulae of Ref. which include complete quark 
mass corrections up to order aj. to evaluate the 
perturbative prediction of R{s) entering into the 
integrals (|) and @). 

Theoretical uncertainties 

Details about the parameter errors used to es- 
timate the uncertainties accompanying the theo- 
retical analysis are given in Refs. [ p4|j29| ]. The- 
oretical uncertainties arise from essentially three 
sources 

{i) The perturbative prediction. The estima- 
tion of theoretical errors of the perturba- 
tive series is strongly linked to its trunca- 
tion at finite order in as . This introduces a 
non-vanishing dependence on the choice of 
the renormalization scheme and the renor- 
malization scale. Furthermore, one has to 
worry whether the missing four-loop order 
contribution ds gives rise to large correc- 
tions to the perturbative series. An addi- 
tional uncertainty stems from the ambigu- 
ity between the results on R obtained us- 
ing contour-improved fixed-order perturba- 
tion theory (FOPTci) and FOPT only (see 
Ref. [|0|). The value as(M|) = 0.1201 ± 
0.0020 is taken, as the average between the 
results from r hadronic decays and the Z 
hadronic width which have essentially un- 
correlated uncertainties [BOl . 



(ii) The quark mass correction. Since a theoret- 
ical evaluation of the integrals (||) and @) 
is only applied far from quark production 
thresholds, quark mass corrections and the 
corresponding errors are small. 

{Hi) The nonperturbative contribution. In order 
to detach the measurement from theoreti- 
cal constraints on the nonperturbative pa- 
rameters of the OPE, the dominant dimen- 
sion D = 4, 6, 8 terms are determined ex- 
perimentally by means of a simultaneous fit 
of weighted integrals over the inclusive low 
energy e+e~ cross section, so-called spec- 
tral moments, to the theoretical prediction 
obtained from Eq. (p^). Small uncertain- 
ties are introduced from possible deviations 
from the PCAC relations (|l|). 

The spectral moment fit of the nonperturba- 
tive operators results in a very small contribu- 
tion from the OPE power terms to the lowest 
moment at the scale of 1.8 GeV (repeated and 
confirmed at 2.1 GeV), as expected from the 
T analyses The value of {^GG) = 

(0.037± 0.019) GeV* found for the gluon conden- 
sate is compatible with other evaluations [ p2| , ^ . 
The analysis proves that global duality holds at 
1.8 GeV and nonperturbative effects contribute 
only negligibly, so that above this energy pertur- 
bative QCD can replace the rather imprecise data 
in the dispersion integrals (^) and (Q). 

Results 

The R{s) measurements and the correspond- 
ing theoretical prediction are shown in Fig. ^. 
The wide shaded bands indicate the regions where 
data are used instead of theory to evaluate the 
dispersion integrals, namely below 1.8 GeV and at 
cc threshold energies. Good agreement between 
data and QCD is found above 8 GeV, while at 
lower energies systematic deviations are observed. 
The R measurements in this region are essentially 



provided by the 772 |34| and MARK I g] Collab- 
orations. MARK I data above 5 GeV lie system- 
atically above the measurements of the Crystal 
Ball 1^ and MDl |^ Collaborations as well as 
above the QCD prediction. 
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Figure 2. Inclusive hadronic cross section ratio in e+e^ annihilation versus the cm. energy -y/s. 
Additionally shown is the QCD prediction of the continuum contribution as explained in the text. The 
shaded areas indicate regions were experimental data are used for the evaluation of Aahad(Afz) ^^'^ '^^'^'^ 
in addition to the measured narrow resonance parameters. The exclusive e^e~ cross section measurements 
at low cm. energies are taken from DM1,DM2,M2N,M3N,0LYA,CMD,ND and r data from ALEPH (see 
Ref. 1^ for detailed information). 



The combination of the theoretical and exper- 
imental evaluations of the integrals yields the re- 
sults il 



Aahad(M|) = (277.8 ± 2.2exp ± 1.4 



±thco. 



X 10" 



-,had 



(695.1 ± 7.5exp ± 0.7thco) X IQ-i" (17) 



with a significant improvement by more than a 
factor of two for Aa{M^), and a 20% better ac- 
curacy on a^^'^ compared to the numbers ([lT|). 

The authors of Ref. improved the above 
analysis in the charm region by normalizing ex- 
perimental results in the theoretically not acces- 
sible region (at least locally) so that they match 
perturbative QCD at safe energies below and 
above the occurrence of resonances. 



The so-renormalized data show excellent agree- 
ment among different experiments which sup- 
ports the hypothesis made that experimental sys- 
tematic errors are completely correlated over the 
whole involved energy regime. The result (after 
correcting for the small top quark contribution) 
reads |38 



Aahad (M|) 



= (276.7 ± 1.7) X IQ- 



(18) 



Another, very elegant method based on an an- 
alytical calculation of the unsubtracted disper- 
sion relation, corresponding to the subtracted in- 
tegral (H), was presented in Ref. Only the 
low-energy pole contribution is taken from data, 
while the contribution from higher energies is cal- 
culated analytically using the two-point correla- 
tion function given in Ref. [HI, and the renor- 
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malization group equations for the running quan- 



tities. This leads to the precise result |39 
Aahad(Afl) = (277.2 ± 1.9) x IQ-^ . 



(19) 



Both numbers (p^ ) an d ([l9| ) are in agreement 
with Aa(M|) from Eq. (|7[). 

(Ill) Constraints from QCD sum rules 

It was shown in Ref. jH] that the previous de- 
terminations can be further improved by using 
finite-energy QCD sum rule techniques in order 
to access theoretically energy regions where lo- 
cally perturbative QCD fails. This idea was first 
presented in Ref. In principle, the method 

uses no additional assumption beyond those ap- 
plied in Section ^. However, parts of the disper- 
sion integrals evaluated at low-energy and the cc 
threshold are obtained from values of the Adler 
D-function itself, for which local quark-hadron 
duality is assumed to hold. One therefore must 
perform an evaluation at rather high energies 
(3 GeV for w, d, s quarks and 15 GeV for the c 
quark contribution have been chosen in Ref. p9|]) 
to suppress deviations from local duality due to 
nonperturbative phenomena. 

The idea of the approach is to reduce the data 
contribution to the dispersion integrals by sub- 
tracting analytical functions from the singular in- 
tegration kernels in Eqs. (||) and (|^), and adding 
the subtracted part subsequently by using the- 
ory only. Two approaches have been applied in 
Ref. first, a method based on spectral mo- 

ments is defined by the identity 

So 

F= f dsR{s)[f{s)-pnis)] 



1 f ds 



27ri J s 

|s|=so 



[Pniso) - Puis)] Dudsis) ,(20) 



with P„(s) = J^dtpnit) and /(s) = 



^had 



as well 



a{0fK{s)/{3Tr^s) for F 

as /(s) = a(0)M|/(37rs(s'- M|)) for F = 
Aahad(M|)[2m^, The analytic functions 

Pn{s) approximate the kernel /(s) in order to 
reduce the contribution of the first integral in 
Eq. (pO|) which has a singularity at s = and 



is thus evaluated using experimental data. The 
second integral in Eq. (^ ) can be calculated the- 
oretically in the framework of the OPE. The func- 
tions Pn[s) are chosen in order to reduce the un- 
certainty of the data integral. This approximately 
coincides with a low residual value of the integral, 
i.e., a good approximation of the integration ker- 
nel /(s) by the p„(s) defined as [p9l 



Pn{s) = X^'^M ^ 



1=1 




(21) 



with the form (1 — s/sq) in order to ensure a 
vanishing integrand at the crossing of the posi- 
tive real axis where the validity of the OPE is 
questioned |^^. Polynomials of order s" involve 
leading order nonperturbative contributions of di- 
mension D = 2(n 4- 1). The analysis is therefore 
restricted to the linear n = 1 case only. 

A second approach uses the dispersion relation 
of the Adler Z?-function 



DfiQ^)^Q^ / ds 



(s + g2)2 



(22) 



for space-like = —q^ and the quark flavour /. 
The above integrand approximate the integration 
kernels in Eqs. (||) and (0), so that the modified 
Eq. (|2^) reads 



So 

F.JdsR^^-is) 



fis)- 



AfQ 



2 1 



4m?. 



Af lDuds{Q^)-Q^Jds 



(S + Q2)2_ 



2\2 



(?3) 



with a normalization constant to be optimized 
for both Aa(M|) and aj^^^. 

In both approaches, a compromise must be 
obtained between uncertainties of experimental 
and theoretical origins. As the subtracted contri- 
bution increases, the experimental error dimin- 
ishes as expected. However, this improvement is 
spoiled by a correspondingly larger theoretical er- 
ror. The procedure followed is designed to mini- 
mize the total uncertainty |29] . 
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Results 

A fit taking into account the experimental 
and theoretical correlations between the polyno- 
mial moments yields for the first (spectral mo- 
ment) approach (hadronic contribution from 2m^ 
to 1.8 GeV) [H 

Aahad(M|)[o,i.8] = (56.53 ± 0.73exp ± 0.39th) x 10"* 



had 
f^M, [0,1.8] 



(634.3 ± 5.6exp ± 2.1th) x 10" 



while the dispersion relation approach gives 
{y/Q^ = 3 GeV) H 

Aahad(M|)[o,i.8] = (56.36 ± 0.70exp ± 0.18th) x 10~* 



,[0,1.8] 



= (632.5 ± 6.2exp ± 1.6th) x 10" 



Only the most precise of the above numbers 
are used for the final results. The above 
theory-improved results can be compared to 
the corresponding pure experimental values, 
Aahad(M|)[o,i.8] = (56.77 ± 1.06) x lO"'^ and 
«Mo,i.8] = (635.1 ± 7.4) X IQ-^", showing clear 
improvement. 

For the charm threshold region only the ap- 
proach ( p3| ) is used giving for the hadronic con- 
tributions from 3.7 GeV to 5 GeV (a 
15 GeV) pi 



Aahad(M|)[3.7 



had 
"m, [3.7,5] 



(24.75 ± 0.84exp ± 0.50th) x 10" 
(14.31 ± 0.50exp ± 0.21th) X 10" 



for which compared with the pure data results, 
Aahad(M|)[3.7,5] = (25.04 ± 1.21) x lO"* and 
(14.44 ± 0.62) X 10^1", only a slight 



^had 

V, [3.7,5] 



improvement is observed. 

5. FINAL RESULTS 

Table ^ shows the experimental and theoreti- 
cal evaluations of Aahad(Af|), a^^'^ and a^^'^ for 
the respective energy regimefjj. Experimental er- 
rors between different lines are assumed to be un- 
correlated, whereas theoretical errors, but those 
from cc and bb thresholds which are quark mass 
dominated, are added linearily. 

According to Table ^, the combination of the 
theoretical and experimental evaluations of the 



^ The evaluation of a^'^'^ follows the same procedure as 

^had 



integrals (||) and (0) yields the final results 



Aahad (M|) 
a-\Mi) 



(276.3 ±l.lcxp± 1.1th) X 10^4 ^ 
128.933 ± 0.015cxp ± 0.015th , (24) 



^had 



(692.4 ±5.6oxp± 2.6th) x 10" 



-10 



a—- = (11 659 159.6 ± 5.6cxp ± 3.7th) x 10" 



-10 



and 



^had 



(187.5 ± 1.7cxp ± 0.7th) X 10 



-14 



for the leading order hadronic contribution to a^. 
The improvement for a{M^) compared to the pre- 
vious results ([TtI ) amounts to 40% and 
about 17% more precise than (17). 



al'"^ IS 



-,SM 

tribution from non-leading order hadronic vac- 
uum polarization summarized in Refs. ^,|l| to 
be a^'^'^iia/Trf] = (-10.0 ± 0.6) x 10-^°. Also 
the light-by-light scattering (LBLS) contribution 
has recently been reevaluated in Refs. and 
H] of which the average (a^^^iLBLS]) = (-8.5 ± 
2.5) X 10^^° is used here. 

Figure ^ shows a compilation of published re- 
sults for the hadronic contributions to a(M|) and 
a^. Some authors give the contribution for the 
five light quarks only and add the top quark part 
separately. This has been corrected for in Fig. 

6. THE MASS OF THE STANDARD 
MODEL HIGGS BOSON 

The new precise result ( p4| ) for a{s) is exploited 
to repeat the global electroweak fit in order to ad- 
just the mass of the Standard Model Higgs boson, 
Mh- The corresponding uncertainty on sin^^vv 
now reaches 5 x lO^'', well below the experimental 
accuracy on this quantity and on nit (see Table 0). 
The prediction of the Standard Model is obtained 
from the ZFITTER electroweak library H and 
the experimental input is taken from the latest 
review The standard value ||] of a{s) yields 



M 



H 



iSAttl) GeV/c^ 



(25) 



while the improved determination (|2J) provides 
log{MH) = 2.02tEJ;23 , 

(26) 



Mh = {lQ5tll) GeV/c' 
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Table 2 

Contributions to AQ!had(Af|), aji^'^ and to a}^^'^ from the different energy regions. The subscripts in the 
first column give the quark flavours involved in the calculation. 



Energy (GeV) 


Aahad(Afl) X 10* 




^had ^ ;^Ql4 


(2m^ - 1.8)„ds 


56.36 ± 0.70oxp ± 0.18th 


634.3 ±5.6oxp± 2.1th 


173.67 ±1.7cxp± 0.6th 


(1.8 - 3.700)„d. 


24.53 ± 0.28th 


33.87 ± 0.46th 


8.13 ± 0.11th 


i/;(lS,2S,3770)c + (3.7-5)„d«c 


24.75 ± 0.84<3xp ± 0.50th 


14.31 ± 0.50exp ± 0.21th 


3.41 ±0.12,xp± 0.05th 


(5 - 9.3)udsc 


34.95 ± 0.29th 


6.87 ± 0.11th 


1.62 ± 0.03th 


(9.3 ~ 12)^tdsc& 


15.70 ± 0.28th 


1.21 ± 0.05th 


0.28 ± 0.02th 


(12 - Oo)udach 


120.68 ± 0.25th 


1.80 ± 0.01th 


0.42 ± 0.01th 


(2mt - oo)t 


-0.69 ± 0.06th 


f» 


f=s 


{2mn - 00)udscbt 


276.3 ±l.lexp± 1.1th 


692.4 ± 5.6exp ± 2.6th 


187.5 ± 1.7cxp ± 0.7th 



Lynn, Penso, Verzegnassi, '87 
Eidelman, Jegerlehner '95 
Burkhardt. Pietrzyk '95 
Martin, Zeppenfeld '95 
Swartz '96 

Alemany, Davier, Hocker '97 
Davier, Hocker '97 
Kiihn, Steinhauser '98 
Groote et al. '98 
Erler '98 

Davier, Hocker '98 



285 290 295 



Barkov et al. '85 
Kinoshita. Nizic, Okamoto '85 
Casas, Lopez, Yndurain '85 
Eidelman, Jegerlehner '95 
Brown, Worstell '96 
Alemany, Davier, Hocker '97 
Davier, Hocker '97 
Davier, Hocker '98 



ixlO-'°) 



Figure 3. Comparison of AahadiM^) and aj^'*'^ evaluations. The values for Aahad(Af|) are taken from 
Refs. ||,|,|5|-0jll|2|,||,|§,|§, while those for ajj'^'* are from Refs. . 



Considering the direct Higgs search currently 
conducted at CERN by the four LEP experi- 
ments, it is worth remarking that the 21 GeV/c^ 
shift between the two indirect determinations is 
almost entirely due to the region 2.0 - 3.7 GcV 
in the evaluation of Aahad(s) where the avail- 
able e^e~ data are systematically higher than 
the QCD prediction (see Fig. ^). In this respect, 
the preliminary results from BES , which are 
more precise than the earlier measurements in 
this energy range, are observed to nicely agree 
with QCD. We are looking forward to more com- 
plete and more precise results in this crucial en- 
ergy region. 



7. SUMMARIZING THE PROCEDURE 
AND ITS JUSTIFICATION 

We have described a three-step procedure to 
improve the evaluation of hadronic vacuum polar- 
isation occuring in the anomalous magnetic mo- 
ments of the leptons and the running of a. By 
far the most rewarding step was to replace poor 
experimental data on e^e~ annihilation cross sec- 
tions in the 1.8 - 3.7 GeV range and above 5 GeV 
by a precise QCD prediction. The justification 
for believing this prediction at the level quoted 
1%) follows from direct tests using experimen- 
tal data. 
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The two assumptions needed for applying QCD 
to this problem are: (i) R{s) can be approximated 
by Rqcd{s) in an average sense only since an in- 
tegral is computed {global quark-hadron duality) 
and (a) perturbative QCD can be reliably used 
at energies as low as 1.8 GeV. 

These hypotheses have been thoroughly tested 
in the study of hadronic t decays [|o[|l) for the 
dominant isovector amplitude, integrating from 
threshold to 1.8 GeV. Using the precise mea- 
surement of Rry+A and moments of the corre- 
sponding spectral function, the nonperturbative 
contributions were found to be smaller than 1%, 
thus enabling to validate the perturbative QCD 
prediction from 1.8 down to 1.0 GeV. Over this 
range, the precision of the test reaches '-^ 1%, go- 
ing down to ^ 2% near 1 GeV. The consistency 
of the QCD description can be expressed through 
the values of as(M|) found in the different cases 
studied: 0.1202±0.0008e:rp±0.0024t?,±0.0010e«oi 
for {V + AJ = 1), 0.1197±0.0018for (V,! ^J), 
0.1207 ± 0.0017 for (A, / = 1) in r decays [|o|, 
and 0.1205 ± 0.0053 for (1/ , / 0, 1) in e+e" an- 
nihilation 

In essence, the nice properties (quark-hadron 
duality and validity of perturbative QCD calcu- 
lations) observed in an a priori critical energy 
region are applied at higher and safer energies, 
where the achieved precision should be 1% or bet- 
ter. 

8. CONCLUSIONS 

This note summarizes the recent effort that 
has been undertaken in order to ameliorate the 
theoretical predictions for a(M|) and a^^'^, cru- 
cially necessary to maintain the sensitivity of the 
diverse experimental improvements on the Stan- 
dard Model Higgs mass, on the one hand, and 
tests of the electroweak theory on the other hand. 
The new value of a^H^^I) = 128.933 ± 0.021 
for the running fine structure constant is now 
sufficiently accurate so that its precision is no 
longer a limitation in the global Standard Model 
fit. On the contrary, more effort is needed to fur- 
ther improve the precision of the hadronic contri- 
bution to the anomalous magnetic moment of the 
muon below the intended experimental accuracy 



of BNL-E821 g, which is about 4 x 10'^° . For- 
tunately, new low energy data are expected in the 
near future from r decays (CLEO, OPAL, DEL- 
PHI, BaBar) and from e+e" annihilation (BES, 
CMD II, DA$NE). Additional support might 
come from the theoretical side using chiral per- 
turbation theory to access the low energy inverse 
moment sum rules ^ and (0). One could, e.g., 
apply a similar procedure as the one which was 
used in Ref. ||5^] to determine the constant Liq of 
the chiral lagrangian. 
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